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For quantum transport through mesoscopic system, a quantum master equation approach is 
developed in terms of compact expressions for the transport current and the reduced density 
matrix of the system. The present work is an extension of Gurvitz's approach for quantum 
transport and quantum measurement, namely, to finite temperature and arbitrary bias voltage. 
Our derivation starts from a second-order cummulant expansion of the tunneling Ifamiltonian, then 
follows conditional average over the electrode reservoir states. As a consequence, in the usual weak 
tunneling regime, the established formalism is applicable for a wide range of transport problems. 
The validity of the formalism and its convenience in application are well illustrated by a number of 
examples. 
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I. INTRODUCTION 

Quantum transport through mesoscopic nanostruc- 
tures has revealed many impressive features associated 
with a number of unique natures such as the quantum 
interferences, discrete levels, and many-body correlations 
0. Depending on the specific systems/problems under 
study, theoretical formalisms have been developed such 
as the Landauer-Biittiker theory and the non-equilibrium 
Green's function (nGF) approach HI, However, gener- 
ally speaking, neither of them implies universal simplic- 
ity in practice, for instance, in treating mesosopic trans- 
port in the presence of many-electron Coulomb interac- 
tion and inelastic scattering with phonons. In particular, 
it is even more difficult to describe the transient processes 
(i.e. time-dependent transport phenomena). 

In some particular cases, a relatively simpler method 
being able to address these issues is the rate equation 
approach 0, 0, IE IE ■ Originally, the "classical" rate 
equation is in certain sense of phenomenological form Q ■ 
Later efforts include its derivation and quantum general- 
ization in the context of the resonant tunneling system, 
based on the nGF quantum kinetic theory as well 
as its modification to describe quantum coherence which 
typically exists in mesoscopic systems 01 ■ partic- 
ular, a microscopic derivation starting with the many- 
particle wavefunction has been presented by Gurvitz et al 
Q. However, an obvious drawback of Gurvitz's approach 
is its limited validity conditions (i.e. the large bias volt- 
age and zero temperature), which largely restricts the ap- 
plicability. Also, they were unable to derive a general for- 
mula in a "system-Hamiltonian-free" form, which means 
the inconvenience that one has to proceed derivation from 
the very beginning for every specific system in practice. 
In this work, we extend Gurvitz's approach to finite tem- 



perature and arbitrary bias voltage, as done in our recent 
work on quantum measurement 8, 9, 10, 11, 12, 13J. In 
particular, we will establish compact expressions for the 
transport current together with the reduced density ma- 
trix, which can serve as a convenient starting point to 
study a variety of mesoscopic transport problems. 

The remainder of the paper is organized as follows. In 
Sec. II, starting with the second-order cummulant expan- 
sion of the tunneling Hamiltonian, formal expressions for 
the transport current and the associated master equation 
are derived. Sec. HI is devoted to a number of examples 
to illustrate the application of the established formalism. 
Finally, in Sec. IV concluding remarks on the approxima- 
tions adopted and the connection with nGF approach are 
presented. In the Appendix, refinement on the cummu- 
lant second-order approximation is self-consistently made 
by including the level broadening effect. 



II. FORMALISM 

Consider the transport setup schematically shown in 
Fig. 1 which is described by the following Hamiltonian 



H 



Hs{al,a^) 



QL—L.B, fik 



H.C.). 



(1) 



a—L,R fik 



Hs is the (mesoscopic) system Hamiltonian, which can 
be rather general (e.g. including many-body interac- 
tion), (a^) is the creation (annihilation) operator of 
electron in state labelled by "/j," , which labels both the 
multi-orbital and distinct spin states of the system. The 
second term is the Hamiltonian of the two electrodes. 
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FIG. 1: Schematic setup for electrical transport through a 
multi-level mesoscopic system. 



The trace in Eq. Q is over all the electrode degrees of 
freedom, leading thus to the equation of motion of the un- 
conditional reduced density matrix of the system. To de- 
scribe the transport problem, we should keep track of the 
record of electron numbers arrived at the collector (emit- 
ted from the emitter and passed through the mesoscopic 
system in between the two electrodes). We therefore clas- 
sify the Hilbert space of the electrodes as follows. First, 
we define the subspace in the absence of electron arrived 
at the collector as which is spanned by the prod- 

uct of all many-particle states of the two isolated reser- 
voirs, formally denoted as i3*^°^ = span{|^'L) (g) 
Then, we introduce the Hilbert subspace "i?^")" ( n — 
1, 2, • • • ), corresponding to n-electrons arrived at the col- 
lector. The entire Hilbert space of the two electrodes is 



which are also termed as emitter (left electrode) and col- 
lector (right electrode), in some places of this work as 
usual. The third term describes tunneling between the 
electrodes and the system. In this paper the electrode 
reservoir electrons are also attached with the index 
to characterize their possible correspondence to the sys- 
tem states. For instance, this will be the typical situation 
in the spin-dependent transport. 

Introducing the reservoir operators Ffj_ — 

J2aJ2k^a^Lk4afik = f + fB.p., WC re-CXprCSS 

tunneling Hamiltonian as 



the 



E 



H.c. 



(2) 



Then, regarding this tunneling Hamiltonian as perturba- 
tion, the second-order cummulant expansion leads us to 
a formal equation for the reduced density matrix 

p{t) = -^^pit) - f dT{£'{t)g{t,T)£'{T)gHt,T))p{t).{3) 

Jo 

Here the Liouvillian superoperators are defined as 
£(...) EE [Hs, {■■■)], ^ [H', {■■■)], and 

g{t,T){---) = G(t,T)(---)G''{t,T) with G{t,T) the usual 
propagator (Green's function) associated with the sys- 
tem Hamiltonian Hs- The reduced density matrix p(t) = 
TrB[pT(t)], and ((• • • )) = TrB[(- • • )pb] with ps the den- 
sity matrix of the electron reservoirs. Notice that Eq. Q 
is nothing but an alternative form of the quantum mas- 
ter equation under the second-order Born approximation. 
The underlining assumption is that the tunneling Hamil- 
tonian is not strong enough, which makes the second- 
order cummulant expansion reasonable. It is known that 
this approximation applies well to most dissipative sys- 
tems in quantum optics. Noticeably, for most transport 
systems, weak tunneling is also the typical regime where 
various forms of golden-rule type theories are adopted. 
In strong tunneling regime, special technique is required, 
which goes beyond the present second-order Born ap- 
proximation, and other golden-rule type theories jl5| . 



With the above classification of the reservoir states, 
the average over states in the entire Hilbert space "B" 
in Eq. Q is replaced with states in the subspace 
leading to a conditional master equation 



p(")(t) 



dTTrB(„, [C'{t)g{t, t) 



xC'ir)g\t,r)pTit)]. 



(4) 



Here p^"'\t) = TTg(,^)[pT{t)], which is the reduced den- 
sity matrix of the system conditioned by the number of 
electrons arrived at the collector until time t. Now we 
transform the Liouvillian operator product in Eq. into 
the conventional Hilbert form: 



c'{t)g{t,T)c'{T)gHt,r)pT{t) 

[H'{t)G{t,T)H'{T)GHt,T)pT{t) 

-G{t, T)H'{T)G^t, T)pT{t)H'{t)] + H.c. 

[/ - //] + H.c. (5) 



To proceed, two physical considerations are further im- 
plemented as follows: (i) Instead of the conventional 
Born approximation for the entire density matrix pT{i) — 
pit) ® Pb, we propose the ansatz pT{t) — X^n /^''"H^) ® 
p^g\ where p^'-* is the density operator of the electron 
reservoirs associated with n-electrons arrived at the col- 
lector. With this ansatz for the density operator, tracing 
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over the subspace "_B("'" yields 



integrated out, leading to 



Tr 



[I] 



(6a) 



Tr 



[II] = 



x[a^G(t,r)4Gt(t,r)p(")] 
x[atG(i,r)a.Gt(t,r)p(")]} 

x[G{t,T)a,G\t,T)p^"^al] 

x[G(t,T)4Gt(t,r)p(")a^] 

+TrB[/I.(r)pi"-'Vi^MW] 
x[G(t,rKGt(t,T)p(«-i)at] 

+TrB[/fi.(r)p^"+^Vl^(i)] 
x[G(t,r)atGt(t,r)p("+i)a^]}. (6b) 



Here we have utilized the orthogonality between states 
in diflFerent subspaces, which in fact leads to the term 
selection from the entire density operator px- (ii) Due 
to the closed nature of the transport circuit, the extra 
electrons arrived at the collector (right reservoir) will 
flow back into the emitter (left reservoir) via the ex- 
ternal circuit. Also, the rapid relaxation processes in 
the reservoirs will quickly bring the reservoirs to the lo- 
cal thermal equilibrium state determined by the chem- 
ical potentials. As a consequence, after the procedure 
(i.e. the state selection) done in Eq. JHJl, the electron 
reservoir density matrices p^'' and p^^"'^'' should be re- 
placed by p^"* , i.e., the local thermal equilibrium reservoir 
state, which leads the reservoir correlation functions in 

Eq. © to be, respectively, {f 1^,(1) fauir)) = ci]X{t-T), 
and (/a^(t)/L(T)) = Ci-I{t-T). Here (• • •) stands 
for TrB[(- • • )p^''], with the usual meaning of thermal 
average. Obviously, {Fl{t)F,{T)) = C^V {t - r) ^ 
j:a=LMCilUt - r), and {F^{t)F^{T)) = ci-\t - r) - 

Under the Markovian approximation, the time 
integral in Eq. is replaced by ^ dr. This 
approximation considerably simplifies the re- 
sult. For instance, substituting the first term 
of Eq. jnj into the time integral of Eq. , we 



have jT^dTd^^\t 
a^[Clt\-C)ai]p^^^){t). 



-iC(t~ 



)at]p(»)(i) 



T)a^[ 

Other terms can be similarly 



-4;V^"-^^a|,-«y"+^)<^]+H.c.}. (7) 



Here 

are 

form of the 

Gayj|/(i£) — 



=L,R. ^"M ■ 

defined 



The 



Ga^i^(±£)a,,, and A 



(±) 



m 

reservoir 



spectral functions Ca^ll{±C) 
terms of the Fourier trans- 
correlation functions, i.e., 
j^^dt&^Mt)e^"^\ We would like 
to remark here that this time integral leads to "exact" 
energy conservation for electron transfer between the 
electrodes and the central system. This "conservation 
law" would cause errors in the near-resonance bias. For 
instance, at zero temperature and for electrode Fermi 
level(s) lower than but very close to certain system level, 
the present "energy conservation law" does not permit 
any electron occupation on the concerned system level. 
Nevertheless, the nGF-based quantum kinetic theory 
allows occupation under the same condition The 
underlying reason is the neglect of level broadening in 
present treatment, whose inclusion is referred to the 
Appendix. 

The "n" -dependence of Eq. is analogical to the 
usual rate equation, despite its formal matrix/operator 
feature. Each term of Eq. Q can be similarly interpreted 
as for the conventional "c-number" rate equation. Com- 
pared with the Bloch equation derived by Gurvitz et al 
0, in Eq. lO is also coupled to which is 

absence from Ref. [3- This feature originates from the 
general nature that Eq. ((T)) is established under non-zero 
temperature and arbitrary (not necessarily large) bias 
voltage. 

With the knowledge of p^^H^), 

one is readily able to 
compute the various transport properties, such as the 
transport current and noise spectrum [l^ . Remarkably, 
we can derive a compact expression for the current which 
is only related to the unconditional density matrix p{t) — 
J2n P^^H^)- The derivation is started with the physical 
observation that the current can be determined by the 
probability distribution function P{n,t) = Tr[p^")(t)], in 
terms of I{t) = edN{t)/dt, where N{t) = I]„«^(":*)- 
Properly treating the summation over "n" and making 
use of the cyclic property under trace, minor algebra 
based on Eq. |[7|) straightforwardly leads to 



I{t) = eY,nTT[p^^Ht) 

71 



AKat)p(t)-fH.c. .(8) 



Here the unconditional density matrix p = p^"^ satis- 
fies an even simpler equation, which can be easily derived 
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by summing up Eq. ^ over "r 



P 



(9) 



Eqs. © and © together with Q constitute the princi- 
pal result of this work, which can serve as a convenient 
starting point to compute transport current under wide 
range of conditions, such as in the presence of many-body 
Coulomb interaction, at finite temperatures and for arbi- 
trary voltages. Moreover, the current expression and the 
associated master equation are free from state represen- 
tation and the specific system Hamiltonian, which there- 
fore holds the merit of unification in its applications. For 
instance, for quantum transport through an interacting 
system, which is usually a challenging problem, one can 
first diagonalize the isolated system Hamiltonian, then 
do the Liouvillian operation easily in the eigen-state rep- 
resentation. In the following, as application of this ap- 
proach we only illustrate a number of simple examples, 
and remain the systematic applications to more interest- 
ing problems to be the subject of forthcoming works. 



III. ILLUSTRATIVE APPLICATIONS 



A. Single Level System 

As a preliminary application of Eq. (|7|), let us consider 
the resonant transport through a single level system. 
Under wide-band approximation for the electrodes, the 
reservoir electron correlation functions read ci*'' {t—r) = 



(efe), where nQ^^(efe) 



= na{ek) is 
(efe) = 1- 

n„(efe). Then the spectral function can be easily carried 
out as 



the Fermi distribution function, and ni 



4±) ^C'W(±£)a = r„ni±)(i?o)c 



(10) 



Here, Ta = 2nga\ta\'^, with ga the density of states of 
the "a" electrode. In the special case of zero temper- 
ature and large bias voltage fj.^ ^ Eq ^ plr, which is 
in fact the applicable condition of Gurvitz's approach 0, 
we simply have A^^-* = r^a, A^"-* = 0, A^^^ r^a, and 
^t^^ = 0. Substituting these into Eq. (TJl yields 



-rLaV^"^a - rflV"^'^a^] + H.c.} . (11) 

To obtain the matrix element form of this equation, let 
us choose the empty (level) state |0) and the occupied 
state |1) as representation basis. Straightforwardly, by 
computing the matrix elements of the terms of Eq. I|ll|) 
one by one, we obtain 



Poo - -^lPqo 



■(") 
Pii 



RPll 



^ LPoa ■ 



This is the result derived by Gurvitz et al under the limits 
mentioned above ItI. 



B. Multi-Level System 

Now we consider the transport through a multi-level 
system as shown in Fig. 1, under arbitrary voltage and 
at finite temperature. The system Hamiltonian simply 

reads Hs = X]^=i ^/j'^Ji'^a'- Also, let us assume that the 
level separation is much larger than the characteristic 
level widths, i.e., — -E^-i] ^ Tl^Th, which leads to 

Sf_,^C^l{t). This as- 



the correlation function c!^u{t) _ wfj_i, 
sumption neglects the interference effect of electron tun- 
neling through different levels, which is significant only 
in the case \Ef^ — < Tl,Tr- Similar to single level 

system, we have A^^^ = T ii{E^)n^^\E^)a^. For this 
simplified model, the reduced system density matrix is 
the direct product of the individual single level density 
matrix, i.e., p = ^^^iPfj., and the steady-state solution of 
the single level density matrix can be easily obtained as 
=Pm|1)m(1| + (1-P/^)|0)m(0|, where |1)^ (|0)^) stands 
for the occupied (unoccupied) state of the pth level, and 
the occupation probability reads 



nLiE^)TLiE^) +nRiE^)TR{E^) 
rL{E,,)+TR{E^) 



(13) 



(12) 



Substituting the obtained A^^^ and p into the current ex- 
pression Eq. (|SJ), we arrive at an expression for the steady- 
state current as 

This result clearly manifests the typical step-like I- 
V characteristics, where each step corresponds to in- 
volving a new level into the conduction by increasing 
the bias voltage , with the standard resonant current 
eTLrBjiTL + TR). 



C. Non-interacting Coupled Quantum Dots 

In the above multi-level system the quantum coherence 
or nature of quantum superposition of system states is 
not manifested, and the result can be obtained via clas- 
sical rate equations. To reveal more clearly the quantum 
nature of the developed formalism, in this subsection we 
consider transport through system of a coupled quan- 
tum dots 0- In this case, the non-diagonal elements of 
density matrix, which have no classical counterparts, will 
appear in the equations of motion and play essential role. 

The Hamiltonian of the coupled quantum dots reads 
Hs — Eia\ai + E2a\a2 + f2(a|a2 -I- a^^i), where each dot 
contains a single resonant level Ei {E2), and the two dots 
are coupled by 17. In principle, for any system the mas- 
ter equation {Tj) or (O can be expressed and solved in the 
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system eigenstate representation. Here, for the coupled 
quantum dots, we would like to present a more elegant 
method in terms of the language of Bogoliubov transfor- 
mation, to explicitly carry out the supcroperators ^i^''. 
To diagonalize Hs , the standard Bogoliubov transforma- 
tion defines a pair of new electron operators as follows: 
bi = uai + va2, and 62 = ua2 — vai. The desired diag- 
onalized Hamiltonian reads Hs = Eib\bi + E2b\b2- The 
diagonalization condition {E2 — Ei)uv + il{u^ — v'^) = 0, 
together with the normalization condition + — 1, 
uniquely determine the transformation coefficients u and 
V, and the eigen-energies read accordingly, Ei — Eiv? + 
E2v'^ + 2Q.UV, and E2 = Eiv'^ +,£'2"^ - 2iluu. Sim- 
ple algebra leads to £"ai = (— i?i)"M6i — (— £'2)"w62, 
and £"a2 = (-£'i)"u6i -I- {-E2Y'ub2. Notice that in 
the wide-band approximation for the electrode reservoirs, 
&^\±C) = Tan;t\-C). We thus have 



4"^ = 



A 



(±) 



Tl 



(^1)61 



(^2)62 



{Ei)bi + un^^\E2)b2 . (15) 



With this result, the explicit form of the master equa- 
tion can be easily obtained for arbitrary offset of the 
dot levels {Ei and £'2)- To compare with the Bloch 
equations derived by Gurvitz ei consider the spe- 

cial configuration of the two dot levels in resonance, i.e.. 
El — E2 = Eq. For this setup, u — —v — l/V^, and 
-El, 2 = Eq ^ il. Moreover, in the large bias voltage 
limit fiL ^ E2,Ei 3> fiR, we simply have A'}^'' = Tlcli, 
A^~'> = 0, = 0, and = ^^02. Substituting 

them into Eq. ^ , an explicit form of conditional master 
equation is obtained as 



4Cp^"'> ~ i {[p(")rLai4 -I- 4rfla2p(' 
-aIp(")rLai - rfla2P^""'^4] + H.c.} . (16) 



In the electron number representation {|1), |2), |3), |4)}, 
which correspond to, respectively, the states of no elec- 
tron in the two dots, one electron in the left (right) dot, 
and two electrons in each dot. Eg J |16() can be precisely 
recast to the result derived in Ref. 0, where the quantum 
coherence nature beyond the classical rate equation was 
particularly emphasized. 



D. Single Level System in the Presence of 
Charging Effect 



The above examples do not involve many-electron 
Coulomb interaction. In this subsection, we consider the 
simplest example of transport through single level system 
in the presence of Coulomb charging effect. The system 
Hamiltonian reads Hs = Y^^{Eo + ^nfi)nfj_. Here the in- 
dex /i labels the spin up ( "t" ) and spin down ( "|" ) states. 



and fl stands for the opposite spin orientation. The elec- 
tron number operator n^^ = aj,a^, and the Hubbard term 
Un^ni describe the charging effect. Obviously, the reser- 
voir correlation function is diagonal with respect to the 
spin indices, i.e., ci^)/(t) = (5^iyCi^^(i). We thus have 



Ai' 



CtJj±C)a^ = C'tliTi^o + Uni,)]a,. (17) 



Moreover, for either spin-up or spin-down electrons, the 

spectral functions ci^),(±i?) are identical to Eq. Hl()|l . 
For present system, the four basis states can be chosen 



as 1 



100) 



|10)u, |3) = |01)Ti, and |4) 



|ll)|j. Also in the limiting case of zero temperature and 
large bias voltage {^l ^ Eq + U > Eq ^ hb.), inserting 
Eq. H17|) into Eq. Q and carrying out the matrix elements 
associated with the above four basis states, we obtain 



•(") 
Pii 

.(n) 
P22 

•(") 

P33 

.(n) 
P44 



(n) 



^{Tr + T'^)P22 

^(r^ + r^)P33 



22 

(") 



"n (") 



(n) 



(18a) 

rUr'\(i8b) 
rUr'\(i8c) 

(18d) 



where F^ = {2nga\ta\'^)E=Eo+u- Satisfactorily, Ea. l(TH|) 
is nothing but the result obtained in Ref. under the 
same limiting conditions. 



E. Interacting Quantum Dot with Zeeman 
Splitting 

In this subsection, we revisit the model studied in 
Sec. III-D, but slightly modify it by allowing for a fi- 
nite spin splitting, i.e., Hs = J2^^=^.liE^^"'io■^^ + T'^^l.'>^■^i)^ 
where the non-zero Zeeman splitting is characterized by 
El — El = A. The transport properties of this system 
has been studied recently by Thielmann et al Ha. by ap- 
plying the real-time diagrammatic technique [l^. Here 
we show our master equation approach can solve this 
non-trivial model in a more transparent way. 

As done previously, we first carry out the commuta- 
Hs,afj,] = -Wf^ttf,, where = E^Sff, + 



Noting that [Hs,W^,] = 0, 
= (-W^;,)"a^. Accordingly, A^t^^ = 



tor £ap = 
we have C' 

c'!^p,{^W^)a^. In the wide-band approximation and 
assuming an energy-independent coupling strength F^ 
(F/j) with the left (right) electrode, explicit expressions 
for A^aJ are obtained as ^^^^^ — ri/jjn^^(W^)a^. In 
the occupation number representation, i.e., |1) = |00)||, 
|2) = |10>Ti, |3) = |01)tx, and |4) = |ll)tx, the master 
equation Eq. can be easily solved, and via Eq. JHJ the 
current can be computed quite straightforwardly. In the 
following, we explicitly carry out the result in different 
voltage regimes. For the sake of being able to obtain ana- 
lytic result, we focus on the limiting case of zero temper- 
ature. Moreover, without loss of generality, we assume 
that the bias voltage makes the Fermi level of the right 
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electrode be always lower than the quantum dot energy 
levels during transport. Therefore, all at the four 
energies, say, E^, Ei, + U, and Ei + U, are zero. 

Regime (i): > E^ + U, E^ + U, E^, E^ > ^r. In this 
high bias regime, the corresponding Fermi functions are 
ni+\E^) = n[+\E,) = n[+\E, + U) = n[+\E^ + U) ^ 
1, and the master equation Eq. Q reads 

Pii = -ZFipii + r7?p22 + ^RP33, 

P22 = -(r_R + rL)p22 + Tipii + ri^p44, 

P33 = -(r_R + rL)p33 + TlPii + rj^p44, 

P4A = -2r_RP44 + rLP22 + Tipas. (19) 

To evaluate the stationary current, only stationary so- 
lution is required, which are easily obtained as, respec- 
tively, P22 = P33 = rLr_R/(rL-i-rfl)2, pn = (r_R/rL)p22, 

and P44 = (ri/ri{)p22. Then, from Eq. (jSJ the current 
is simply carried out as 



lit 



e^R{p22 + P33 + 2p44) 



2erLr/ 



(20) 



Regime (ii): E-^ + U > fiL > E^ + U,E^,Ei > pij. 
The Fermi functions in this case read n^^''(£'| + U) = 



0, n[+\E^) = nriEi) 
resulting master equation is 



(+)/ 



n^j^\E, + U) 



1, and the 



Pll = -'i.TLPu +TrP22 +TrP33, 

P22 = -(r_R + rL)p22 + Elph + r7?p44, 

P33 = -^RP33 + TlPii + r7^p44 + rLP44, 

P44 = -(Fl + 2rj^)p44 + rLP22- 



(21) 



Solution of the stationary state reads, respectively, P22 = 
TLTRiTL + 2Tr)/2{Tl + rfi)3, P33 = Pll + 2ri (Fi + 
rR)/2{TL + Tr)^, and p44 = rfl/2(rL + Tr^. Note 
that Pll = 1 — P22 — P33 — P44, which is irrelevant to the 
current. Straightforwardly, we obtain the current 



I{t ^ 00) = er_R(P22 + P33 + 2p44) 

^ er^rfl(r^ + 2rfl) 

{Tl+Tr)^ ■ 



(22) 



Regime (Hi): E^ + U, Ei + U > pi > E-f,Ei > pR. 



The Fermi functions n^L ^^l) = "-r^(^T) = 1' and 



i^l\e^ + U) = n^L'{Ei + U) = 0. The corresponding 
master equation reads 

Pll = -2FlPi1 + F_RP22 + F/JP33, 

P22 = -^RP22 + TlPii + (Fl + Fi^)p44, 

P33 = -^RP33 + TlPii + (Fl + Ffl)p44, 

P44 = -2(Fl+F;j>44. (23) 

The stationary-state solution of the reduced density ma- 
trix leads to the transport current as 



.(+) 



I(t ^ 00) = eFii;(p22 -I- P33 + P44) 



2eFLFj^ 
2FLH-Ffl' 



(24) 



Regime (iv): E^ + U,Ei + U, E^ > pL > E^ > pR. In 

this setup, only n^j^\E^) — 1, and all other Fermi func- 
tions are zero. Similarly, we first carry out the master 
equation 

Pll = -FlPiI + rLP22 + F7?P22 -I- F7^P33, 

P22 = -(Fl + r7^)p22 + (Fi + Fi?)p44, 

P33 = -ri?p33 + Fipii + {Tl + Tr)p44, 

P44 = -2(Fl+Fk)P44. (25) 

Then, the stationary transport current is calculated via 
the stationary-state solution of the density matrix as 



I{t ^ 00) = eF_R,(p22 -I- P33) = 



eFiFfl 



(26) 



Remarkably, we have precisely recovered all the 
Coulomb plateaus presented in Ref . 0, which go beyond 
simple intuition and are obtained there by a not easily 
accessible real-time diagrammatic technique. This exam- 
ple may shine light on the convenience of our approach 
in applications. 



IV. CONCLUDING REMARKS 

In summary, we have developed an efficient master 
equation approach for quantum transport through meso- 
scopic systems, and demonstrated its application by a 
number of examples. Compared with the previous work 
by Gurvitz 7], the present study not only generalizes 
the applicable conditions to finite temperature and ar- 
bitrary voltage, but also identifies the adopted approx- 
imation which appears not very clear in Ref. 0. That 
is, by treating the electrodes as (Fermi) thermal baths, 
the major approximation adopted in our derivation is 
the standard second-order Born approximation for the 
coupling (tunneling) Hamiltonian. It is known that this 
well-justified approximation makes the resultant quan- 
tum master equation applicable in a large number of dis- 
sipative systems (e.g. in quantum optics), provided the 
system-bath coupling is not so strong. Favorably, the il- 
lustrated examples in this paper also show its applicabil- 
ity in quantum transport. Moreover, the developed mas- 
ter equation approach holds the obvious advantages of 
application convenience and straightforwardness, as well 
as the ability to address many-body correlation, inelas- 
tic scattering, and transient behavior, which are usually 
difficult issues in mesoscopic transport. 

In comparison with the nGF approach, we found that 
the structure of Eq. |(SJ is in fact identical to the formal 
expression of current in terms of non-equilibrium corre- 
lation functions ^] . The nGF approach remains the rel- 
atively hard task to searching for particular techniques 
(e.g. the Feynman diagram or equation of motion) to 
carry out those correlation functions. In this sense, the 
obtained Eq. ijHJl is nothing but the explicit Markovian re- 
sult under the second-order Born approximation for the 
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tunneling Hamiltonian. In principle, further systematic 
corrections are possible along the line of going beyond 
the Born approximation, to include higher order contri- 
bution of tunneling. Finally, we mention that Eq. ^ can 
be derived from the formal nGF expression of current, 
however, the present derivation along the line of Ref. Q 
is interesting, and the particular result Eq. jT)) from this 
unique method is of great value, which contains rich infor- 
mation and can be conveniently employed, for instance, 
to calculate the noise spectrum 0, 0, 0| . 

APPENDIX: LEVEL-BROADENING EFFECT 

To make the derived formulas Eqs. Q)-® more accu- 
rately applicable to arbitrary voltage, additional care is 
needed as any of the individual system levels is approach- 
ing to the Fermi surface of the electrode. For the sake 
of description clarity, let us take the single-level reso- 
nance system as an example to highlight the key point. 
As mentioned previously, the treatment in Sec. II un- 
der the second-order Born approximation has neglected 
level broadening effect, which would cause certain errors 
in some particular cases. For instance, current flowing 
through the resonance system would be strictly forbid- 
den under the near-resonance condition, i.e., as the res- 
onance level Eq is a little bit higher than /i^ (/i/, > ^n). 
However, it is well known that a full quantum treatment 
will give nonzero current in this situation 0, Q . For bias 
voltage such that the resonance level (together with its 
broadening) is within the range of the two Fermi levels 
(i.e. fj,L > Eq > ^j-r), common result of resonance cur- 
rent would be predicted by our master equation approach 
and the nGF-based quantum transport theory. In spite 
of this, it would be desirable to remove the drawback of 
inaccuracy of our approach in the near-resonance situa- 
tion. 

To account for the level-broadening effect, we return 
to the evaluation of a'^^ — ci^t(±£)ay . Without 
loss of generality, we restrict our description in the diago- 
nal case C^i(t) — S^i,C^^fj,{t). More general description 
is straightforward provided one has clarified the correla- 
tion between "/^" and "/i^" . Using the free-electron-gas 

model for the electrodes, Aq^-* can be expressed as 
Aii^ = 2^ |i„^,,pni±)(6fc) / die±-^-*e±'^*a^.(A.l) 

In our previous treatment, we have replaced the time 
integral 2 dt by dt, under the spirit of Marko- 



vian approximation. As a result, the time integration 
gives rise to a ^-function, 2'K5(tk + £), which charac- 
terizes energy conservation for electron transfer between 
the central system and the electrodes. Mathematically, 
this procedure is equivalent to dropping the imaginary 
(principal) part of integral, and keeping only the real 
part. Now, notice that e^*'^*a^ describes the quantum 
evolution of the E^ state (level) associated with the iso- 
lated system Hamiltonian. As a standard procedure, the 
system level broadening effect due to coupling with the 
electrodes can be implemented by inserting a damping 
factor e"^*"* = e"^^^''"'"^"''^* into the integrand of the 
time integral After this, by keeping only the real 

part of the integral and adopting the typical wide-band 
approximation for the electrodes, we have 

A'^aJ = T^m/ '^ao.,{eu+C)nt\eu)a,. (A.2) 

Here the standard Lorentzian spectral density function 
reads aa^(w + C) = 2r^/[(tj + Cf F^]. Formally in- 

troducin g Ni '^\-C) = / |^aQ^(efc + C)n'i^\ek), we re- 
express (jA.2|) in a very compact form as 

^L*^' = raM^i?(-^K- (A.3) 

Elegantly, n!^^ i^^) can be regarded as the counterpart 
of the Fermi function nl^\—C) after accounting for the 
level broadening. Combining (|A.3|1 with Eqs. Q)-©, we 
complete the generalization of the formalism. 

As an illustrative application of the generalized for- 
malism, we re-consider the transport through the (free) 
multi-level system. Straightforwardly, the current ex- 
pression of Eq. H14|) becomes 

^ ^ ^ rL{E^)rR{E^) 

Z^rL{E^)+Tn{E^) 

-^aaf,{ek - E^)[nL{ek) - n^Xefc)]- (A.4) 

This is the well-known formula for resonant tunneling 
current, which is valid for arbitrary voltage including the 
near-resonance situation. 
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